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We report on recent results of point-contact spectroscopy measurements in Mg^^B2 polycrystals 
irradiated at different neutron fluences up to $ = 1.4 • 10^°cm~^. The point contacts were made by 
putting a small drop of Ag paint - acting as the counterelectrode - on the cleaved surface of the 
samples. The gap amplitudes were extracted from the experimental conductance curves, showing 
Andreev-reflection features, through a two-band Blonder- Tinkham-Klapwijk fit and reported as a 
function of the Andreev critical temperature of the junctions, The resulting Ao■(^J^) and 

^tt(T^) curves show a clear merging of the gaps when T^ ~ 9 K that perfectly confirms the 
findings of specific-heat measurements in the same samples. Anomalous contacts with T^ > Tc 
(being Tc the bulk critical temperature) were often obtained, particularly in samples irradiated at 
very high fluences. Their fit gave a different dependence of A,r on T^ . The possible origin of these 
contacts is discussed in terms of local current-induced annealing and/or nanoscale inhomogeneities 
observed by STM in the most irradiated samples. 

PACS numbers: 74.45.-l-c, 74.70.Ad, 74.62.Dh 



I. INTRODUCTION 

Most of the present fundamental research on the two- 
band superconductor MgB2 is devoted to studying the 
effects of substitutions and disorder on its properties. 
This interest in exploring the "neighborhood" of the pure 
compound is justified in part by the quest of a recipe for 
improving some of its properties - especially in view of 
power or electronic applications - and in part by the need 
of understanding at the best this unique example of two- 
band phonon-mediated superconductor with a relatively 
high Tc- As a matter of fact, the presence of two systems 
of bands crossing the Fermi surface, each developing an 
energy gap below T^, has a number of intriguing con- 
sequences that make the physics of MgB2 unexpectedly 
rich and complex. One of these aspects is the role of scat- 
tering by impurities. Due to the different parity of the a 
and TT bands, scattering of quasiparticles between them is 
highly improbable in the pure compound, while almost 
independent scattering rates exist within the bands. In- 
traband scattering has no effect on the gaps or on , but 
affects various properties of the material, e.g. the critical 
fieldJ and the magnetic-field dependence of the gaps^i. On 
this basis, the ratio of the diffusivities in the two bands 
has been experimentally evaluated in pure MgBa^*^. 

According to an early prediction of the two-band 
model&j the increase of interband scattering in a system 
like MgB2 should make the two gaps approach each other 
and finally merge into a single BCS gap. However, ob- 
serving this effect in a real material has turned out to be 
more difficult than expected. The few chemical substitu- 
tions that actually take place (e.g. C substitution for B, 
Al or Mn substitution for Mg) give rise to lattice or elec- 
tronic effects that can mask the increase in disorder. For 



example, C and Al substitutions also cause remarkable 
changes in the DOS at the Fermi level, in the phonon 
frequencies, in the cell volume and so on, with an ob- 
vious complication in the interpretation of the data. In 
C-substituted single crystals the merging of the gaps has 
been recently observed"* as arising from the interplay of 
the band filling due to electron doping and an increasing 
amount of interband scattering, probably due to extrinsic 
reasons 

The controlled damaging of the compound by means 
of irradiation allows partly overcoming these difficulties, 
even though the irradiated material is far from being an 
"ideal" disordered version of MgB2. In particular, side 
effects of irradiation are the creation of Li atoms and He 
nuclei in the lattice (due to the thermal neutron capture 
by ^"B nuclei), the reduction in the partial DOS of the 
2px,y statesifi, the anisotropic expansion of the crystal 
latticeii. The effect of neutron irradiation (up to very 
high fluences) on the energy gaps of MgB2 has been re- 
cently studied by means of specific- heat measurements™, 
showing the achievement of single-gap superconductivity 
in samples with Tc as low as 11 K. 

In this paper, we present and discuss the re- 
sults of point-contact spectroscopy measurements (in 
the Andreev-reflection regime) in the same neutron- 
irradiated Mg^^B2 samples studied in Ref. We will 
show that the gap amplitudes measured by PCS agree 
very well with those given by specific-heat measurements 
and we will discuss the experimental trend of the gaps 
within the two-band Eliashberg theory. 

Moreover, we will report on the anomalous features 
of a large number of contacts whose Andreev critical 
temperature is greater than the bulk Tc. These 
contacts feature very good Andreev-reflection conduc- 
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tance curves that were very well fitted by the two-band 
Blonder-Tinkham-Klapwijk (BTK) model to extract the 
gap amplitudes Ac, and A^r- Once reported as a function 
of T^, A^ has a completely different trend with respect 
to that reported in Ref. [13 and observed by PCS in 
"standard" contacts. We will discuss this odd result in 
terms of local nanoscale inhomogeneities of the material 
and/or local annealing due to the technique we used to 
tune the properties of our "soft" point-contact junctions. 



II. EXPERIMENTAL DETAILS 

The procedure for sample fabrication and irradiation 
is reported in detail elsewher o^^i^^ . The samples were 
prepared by direct synthesis from pure elements, using 
in particular isotopically-enriched ^^B (99.95% purity) 
with a residual ^"B concentration lower than 0.5% so 
as to make the penetration depth of thermal neutrons 
greater than the sample thicknessiS. The samples we 
measured had been irradiated at the Paul Scherrer Insti- 
tute (PSI) in Villigen, Switzerland. For simplicity and 
ease of comparison, let us label them as in Ref, 12, i.e. 
PO (pristine Mg^Bz), P3 (fluence $ = 7.6 • lO^^cm-^), 
P3.7 ($ = 5.5 • IQi^cm-^), P4 ($ = 1.0 ■ lO^^cm-^) and 
P6 ($ = 1.4 • 10^°cm~^). Many of their structural and 
transport properties are reported in Refs. and IT^ 
The bulk critical temperatures, defined as Tc = Tqq% of 
the superconducting transition measured by susceptibil- 
ity, are: 38.8 K for PO, 35.6 K for P3, 25.8 K for P3.7, 
20.7 K for P4, 8.7 K for P6. The width of the supercon- 
ducting transition, defined as Arc(10% — 90%), varies 
from 0.3 K (in PO and P3) up to a maximum of 0.9 K (in 
sample P4)i^. The transition remains rather sharp also 
in the most irradiated sample, which indicates a highly 
homogeneous defect distribution even at the highest flu- 
ence. This homogeneity is also confirmed by the sharp 
X-ray diffraction peaks for the (002) and (110) reflections 
reported in Ref. Il3l which also indicate an anisotropic 
expansion of the cell parameters, more pronounced along 
the c axis (up to 1%). The residual resistivity increases 
by two orders of magnitude (from 1.6 /ifJ-cm for PO to 
130 /if2-cm for P6) with a corresponding reduction in the 
residual resistivity ratio (RRR). 

The point contacts were made by placing a small 
drop of silver paint on the freshly cleaved surface of the 
samplesii. The conductance curves, dl/dV vs. V, were 
obtained by numerical differentiation of the measured 
I — V curve. In all cases, we studied the temperature 
dependence of the curves, which show clear Andreev- 
reflection features, so as to determine the critical temper- 
ature of the junction (in the following referred to as the 
"Andreev critical temperature", T^). Strictly speaking, 
in fact, rather than the bulk Tc is the critical temper- 
ature to be related to the local gap amplitudes measured 
in a given contact. In point-contact spectroscopy, can 
correspond to any temperature between the onset and 
the completion of the superconducting magnetic transi- 



tion. Therefore, one usually has = Tc within the 
experimental broadening of the superconducting transi- 
tion, i.e. Tq% < T^ < rioo% (let us recall that here we 
defined Tc = TgQ%). In irradiated samples, this actually 
occurs in a subset of contacts we will call "standard" 
contacts. The conductance curves were divided by the 
normal-state conductance and then fitted with a two- 
band BTK model in which the conductance through the 
junction is expressed by G = (1 — uj-,r)Ga 4- WttGtt, 
and being the partial a- and 7r-band conductances, 
and bJ-jr the weight of the 7r-band contributiorjiiii^. The 
model contains as adjustable parameters the gap ampli- 
tudes Ao- and A^, the barrier parameters and Z^, 
the phenomenological broadening parameters and F^, 
plus the weight w^^ . The broadening parameters enter in 
the definition of the density of states in the usual way, 
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In this context, they account for both intrinsic (i.e. fi- 
nite lifetime of quasiparticles) and extrinsic (related to 
the technique and the nature of the contacts) phenom- 
ena that smear out the conductance curves. Zcr,TT are re- 
lated to the potential barrier height at the interface and 
to the mismatch in the Fermi velocity vp between the 
two sides of the contact. Owing to the different values 
of v-p in the a and tt bands, we allow Z„ ^ Z^^. Finally, 
Wtt is predicted to range from 0.66 to 0.99 for perfectly 
directional tunneling in pure MgB2^, depending on the 
angle of current injection with respect to the ah planes. 
In the absence of specific predictions in samples with re- 
duced anisotropy, we kept w^^ in the same range as in 
Refs. im 

We generally selected contacts with rather high values 
of the normal-state resistance i?N , corresponding to small 
values of the contact size a which has to be smaller than 
the electronic mean free path I for energy-resolved spec- 
troscopy to be possible. The limit condition a ^ £ de- 
fines the so-called ballistic regime of conductioni^. Based 
on Ref. la we also required the conductance curves of 
our point contacts not to present dips, which are the 
hallmark of a breakdown of the conditions for ballistic 
conduction at finite voltage and signal the presence of 
heating in the contact region. When the contact resis- 
tance was too small, or its conductance did not show 
clear Andreev-reflection features, we were able to change 
the contact characteristics (in a surprisingly repeatable 
way) by applying short voltage or current pulses to the 
junction itself. In some cases, we also used the magnetic 
field to clarify whether one or two gaps were present, as 
explained in detail elsewherei^ and in the following. 



III. RESULTS IN STANDARD CONTACTS 

Fig. ^ reports an example of the raw conductance 
curves measured as a function of the temperature in three 
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point contacts made on samples P3, P4 and P6. The 
low-temperature curves clearly show the typical Andreev- 
reflection features - in particular, the two symmetric 
maxima at ibV^oak approximately corresponding to the 
edges of the small gap. The large gap-features are usually 
less clear even in pure MgB^^ and in disordered samples 
they are difficult to see. The thick curve in each panel 
indicates the normal-state conductance and the relevant 
temperature is thus chosen as the of the contacts. 
Note that, in all cases, lies between the begin and 
the end of the magnetic superconducting transition. 

In Fig.^a) and (b) the shape of the normal-state con- 
ductance curves and the absence of dipsi^ in the super- 
conducting state indicate that no heating occurs in the 
junctions, which are thus in the ballistic regime (i.e., 
the contact size a is smaller than the mean free path 
£) and energy-resolved spectroscopy is possible. Here 
and in the following £ is defined through the relation 
= £~^ + £^^ , where £c and £[ are the elastic and in- 
elastic mean free paths, respectively (at sufficiently low 
temperature, £ ~ £o)- 

In Fig. ^c) small dips at \V\ > |Vpcak| suggest that 
the contact might actually be in the diffusive regime 
(i.e. -^o < a < A, being A the diffusion lengthil.) in 
which energy-resolved spectroscopy is still possible. In 
this regime, the resistance of a point contact between 
two normal metals - in the hypothesis that the Fermi ve- 
locity of the metals are almost equal and the resistivity of 
one metal (here Ag) is much smaller than the resistivity 
of the other (here sample P6) - can be expressed hyi^ 
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where the first term is the Sharvin resistance for ballistic 
conduction and the second one is the Maxwell resistance 
for the thermal regime multiplied by a function of the 
Knudsen ratio k — £/a. This function, F, is always of 
the order of unity, p is the normal-state resistivity of the 
irradiated sample (that we will identify with the resid- 
ual resistivity) and £ is the mean free path (evaluated in 
Rcf 13). In the second term, the contribution of the first 
half of the contact (the Ag counterelectrode) has been 
neglected^fliSi due to the much smaller resistivity of Ag 
with respect to that of sample P6. At T < Tc the irra- 
diated MgB2 is superconducting and, thus, the contribu- 
tion of the second term in eq.^should disappear. Never- 
theless, it is easy to show that even at very low bias volt- 
ages - of the order of the energy gap, here about 1 meV 
as we will discuss later - the current density in the con- 
tact is higher than the critical current density and thus 
tends to drive normal a small volume of the superconduc- 
tor pushing back the NS boundary a short distance^Si. If 
the size of this normal region is smaller than the coher- 
ence length the spectroscopy of the gap is still possible^Si, 
but the second term in eq. ^ starts playing a role and 
a small, voltage-dependent and temperature-dependent 
heating appears in the contact. The small vertical shift 
in the conductance curves shown in Fig. ^c) on increas- 
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FIG. 1: Temperature dependence of the raw conductance 
curves measured in three contacts on samples P3, P4 and 
P6, whose bulk critical temperature is indicated. Thicker 
lines indicate the normal-state conductance, which, in most 
of the contacts, is practically temperature-independent and is 
reached when T = T^^. The temperature of each curve is also 
indicated in the legend. 



ing temperature confirms this picture indicating that the 
temperature-dependent resistivity of the material plays 
a role in the contact resistance. Heating of the contact 
region generally gives rise to an apparent decrease of the 
critical temperature of the contact with respect to 
the bulk Tc. However, here the superconducting features 
disappear at some temperature between 8.0 and 8.5 K, 
which is only slightly smaller than Tc = 8.7 K. Hence, we 
can conclude that a very moderate heating is likely to oc- 
cur in the contact shown in Fig.^c) and it can be safely 
neglected as long as the voltage drop across the junction 
is of the order of Fpcak- We will show in greater detail 
in a following section that the two conditions described 
above (normal-region size < ^ and very small heating) 
are compatible with the curves shown in Fig.^c) only if 
parallel diffusive nanocontacts are supposed to be present 
in the contact region. 

Fig. 121 reports an example of experimental, normalized 
conductance curve (symbols) for each sample. Notice 
that the horizontal scale is the same for all the panels, 
so as to highlight the shrinking of the Andreev-reflection 
structures on increasing the neutron fluence - which in- 
dicates, in turn, a decrease in the amplitude of the gaps. 
While in the top curve (sample PO) the presence of peaks 
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FIG. 2: Normalized conductance curves (symbols) of different 
point contacts on neutron- irradiated Mg^^B2 polycrystals at 
T=4.2 K. The curves are labeled with the name of the samples 
and the relevant bulk Tc according to Ref. 17^ The values of 
the normal-state resistance are also indicated. Solid lines are 
the best-fitting curves given by the appropriate BTK model 
(two-band in samples PO to P4, single-band in sample P6). 
Dashed lines shown for samples P3, P3.7 and P4 represent 
the single-band best-fitting curve to be compared to the two- 
band fit. In the case of sample P4, solid and dashed lines are 
almost superimposed. 



and shoulders clearly witnesses the existence of two gaps, 
in the irradiated samples this evidence is lacking, due to 
a progressive broadening of the curves accompanied by 
a reduction in their height. The same happens in doped 
MgB^LSi. In all these cases, the existence of two gaps can 
be inferred from the fit of the curves with the BTK model 
or evidenced by the application of a magnetic field^. 

The BTK curves that best fit the experimental data 
are shown in Fig. |21 as solid lines. In sample PO (pris- 
tine Mg-'^^B2), the fit can only be obtained with the two- 
band BTK model. In samples P3 and P3.7, the two-band 
fit works better than the single-band one, since it repro- 
duces both the width of the Andreev-reflection structures 
and the position of the peaks, while the single-band fit 
(dashed lines) does not. The curve measured in sam- 
ple P4 admits both a single-band and a two-band BTK 
fit, which are almost equally good - as a matter of fact, 
dashed and solid lines are almost superimposed in this 
case. In sample P6, the dips at \V\ > |lpcak| modify the 
shape of the curve so that asking the model to fit the 
curve in this region is nonsense. In these conditions, the 
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FIG. 3: Black symbols: gap amplitudes from PCS mea- 
surements, as a function of the Andreev critical tempera- 
ture. Open symbols: gap amplitudes from specific-heat 
measurements-- as a function of the critical temperature mea- 
sured by specific heat. Solid lines are only guides to the eye, 
while the straight dashed line indicates the value of the gap 
for a conventional, BCS superconductor. 



parameters of the BTK model should be adjusted so as 
to fit the conductance maxima and the zero-bias dip be- 
tween them. The single-band BTK model is sufficient to 
accomplish this task very well (see the line in the bottom 
panel of Fig.O. If a two-gap fit is tried, the values of 
and Ao- turn out to be so close to each other to be prac- 
tically identical. Hence, in this sample the existence of 
a single gap can be safely concluded, in agreement with 
the findings of Ref. 12. 

The gap amplitudes extracted from the fit of the curves 
shown in Fig. [3 (and of other curves not reported here, 
measured in different contacts on the same samples) are 
plotted in Fig.|3|as a function of (black symbols). In 
the region around 18-19 K, the gap amplitudes result- 
ing from the two-band fit of the conductance curves are 
shown, but it should be borne in mind that a single-gap 
fit is possible as well in this region, giving a gap value 
A ~ A^. In the same figure, the gap amplitudes ex- 
tracted from the fit of specific-heat measurements^^ are 
also shown (open symbols). The agreement between the 
two sets of data is good in the whole range of critical 
temperatures, especially if one takes into account that: 
i) PCS is a local, surface-sensitive technique while spe- 
cific heat is a bulk property; ii) the gap values obtained 
by PCS are correctly plotted versus the Andreev critical 
tem pera ture of the contacts, T^, while those taken from 
Ref. Il3 are reported as a function of the specific- heat Tc- 

The trend shown in Fig.|31clearly indicates a transition 
from two-band to single-band superconductivity at high 
neutron fluences. That the heavily-irradiated material 
undergoes deep changes above a certain neutron fluence 
(<& ~ 10^^ cm~^) is confirmed by the steep decrease in 
Tc and Bc2, by the increase in the cell parameters a and 
(jiiii^ and by the observed decrease in the tr-band DOS^°. 
Incidentally, it is interesting to note the initial, small 
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FIG. 4: (a) Symbols: gap amplitudes from PCS measure- 
ments, as a function of the Andreev critical temperature, T^. 
Lines: gap amplitudes vs. calculated by solving the two- 
band Eliashberg equations using as adjustable parameters: i) 
the interband scattering F^tt alone (dashed lines); ii) the a- 
band DOS A^o-(-Bf) alone (thin solid lines), (b) The values of 
Fo-TT that give rise to the dashed curve of panel (a) (previous 
case (i)). (c) The NaiEp) vs. Tc curve necessary to fit the 
data (previous case (ii)). 



increase in A^r, which is the hallmark of an increase in 
the scattering between bands. 

We tried to reproduce the experimental trend of the 
gaps and reported in Fig. |21 within the two-band 
Eliashberg theory. A complete fit of the ATr{T^) and 
Acr{T^) curves is actually impossible since in a certain 
range of critical temperatures both gaps are smaller than 
the BCS value (see Fig. until a BCS-like gap ratio is 
almost recovered at Tc ~ 9 K, when A ~ 1 meV. An 
energy gap smaller than the BCS value has indeed been 
observed in disordered, conventional superconductors^^ 
and the same might occur in a two-band system, but it 
is strictly forbidden within the Eliashberg theory and no 
explanation for these findings has been given yet. Once 
established this point, one can proceed with the fit. 

The simplest approach is to consider the irradiated ma- 
terial as if it was only "disordered MgB2" , thus neglecting 
the changes in the DOS, in the phonon frequencies and in 
the cell volume, and only increasing the interband scat- 
tering r^^. Once the value of this single parameter is 
chosen to reproduce the critical temperature of a given 
sample, no further degrees of freedom are left to repro- 
duce the gap amplitudes. The resulting curves are shown 
in Fig. ^a) as dashed lines and the corresponding values 
of Fo-TT are reported in Fig. 01b). Notice that very high 
values of Fq-tt would be necessary to suppress Tc below 
30 K and make the gaps merge. Such an intense inter- 
band scattering is never observed in real systems and 
is probably not physical. In doped MgB2, for example, 
the suppression of Tc and A^ is mainly due to other ef- 
fects (typically a reduction in the DOS^) and even if F^tt 
moderately increases (usually remaining smaller than 10 
meV) its effects are partially masked, so that a tendency 
of Att to remain constant or slightly increase is at most 



observed2iS. The present case is not very different since, 
as previously pointed out, neutron irradiation has "side 
effects" such as sizeable changes in the a DOS and in the 
cell parameters that are not included in this description. 

The opposite approach for the fit of the experimental 
gap values of Fig^Ja) thus consists in disregarding the ef- 
fect of disorder (scattering) and only taking into account 
the change in the a DOS at the Fermi level, Na{EF)'^°. A 
reasonable fit of the experimental Acr(Tc'^) and A^(Tc'^) 
curves (with the general theoretical limitation that the 
gap ratios cannot be both smaller than the BCS one) 
is indeed obtained in this way, as indicated by the solid 
lines in Fig. ^a). That using a single parameter one 
can reproduce in such a good way the values of the two 
gaps and the critical temperature is, by itself, a good 
result and indicates that N^{Ep) is largely dominant in 
determining the observed gap trend. This conclusion is 
consistent with the findings of Ref. where the depres- 
sion of Tc down to about 10 K was justified by inserting 
in the McMillan formula the reduced DOS (about 25% of 
the value in pristine MgB2) measured by NMR. Fig.^fc) 
reports the Tc^ dependence of the cr-band DOS necessary 
to fit our PCS data (solid fine). 

To account for the initial increase in A^ (also clearly 
shown by specific-heat measurementsi^) a small amount 
of interband scattering must be inserted in the model. 
In particular, in our case Fo-tt should be about 0.7 meV 
when Tc=35 K and should saturate at a constant value 
(no more than 2 meV) at low Tc. It is clear, however, 
that even including in the model all the possible effects 
of neutron irradiation, the agreement with the data can 
be hardly improved due to the aforementioned anomaly 
of the gap values that are both smaller than the BCS 
value in the Tc^ range between ~ 9 K and ~ 20 K. 



IV. RESULTS IN ANOMALOUS CONTACTS 

The percentage of "standard" contacts (as defined in 
the previous section) is equal to 100% in samples PO and 
P3, but fast decreases in more irradiated samples. For 
example, it is about 70% in sample P3.7, 30% in sam- 
ple P4 and becomes as small as 10% in sample P6. The 
remaining contacts are "anomalous" in the sense that 
their Tc^ exceeds Tc, which is clearly related to some 
kind of intrinsic or induced inhomogeneity in the sam- 
ples. Fig. O reports the conductance curves of one of 
such anomalous contacts on the most irradiated sample 
(P6, bulk Tc = 8.7 K). The normal-state resistance of 
the contact was i?N — 71 fl. The temperature depen- 
dence of its conductance curve, reported in Fig. EJa), 
clearly shows that the Andreev-refiection features persist 
well above the bulk critical temperature and disappear at 
Tc"* = 32.7 K, which is more than three times the bulk Tc 
measured by susceptibility. The thick line in Fig.[Sfa) in- 
dicates the normal-state conductance curve. A fit of the 
low-temperature curve with the BTK model unambigu- 
ously shows the presence of two gaps whose values are 
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FIG. 5: (a) Temperature dependence of the raw conductance 
curves measured in a contact with r^=32.7 K obtained on 
sample P6 (bulk rc=8.7 K). The thick line is the normal-state 
conductance, (b) Magnetic-field dependence of the same con- 
ductance curve as in (a). Again, the thick line is the normal- 
state conductance, (c) Comparison of the conductance curves 
in zero field and in a magnetic field of 2 T. The outward dis- 
placement of the conductance peaks is highlighted by vertical 
lines. 



similar (but not identical) to those obtained in standard 
contacts with the same . 

To further enlighten this point, we applied to the junc- 
tion a magnetic field perpendicular to the direction of 
main current injection, and studied the behavior of the 
conductance curves on increasing the field intensity. In 
pure MgB2 , using a magnetic field allowed us to separate 
the partial contributions of the a and tt bands to the 
conductance across the junctioni2*i^. In doped MgB2, 
the complete separation is not always possible but, if 
two gaps are present, an outward shift of the conduc- 
tance maxima occurs when the smaller gap is strongly 
suppressed by the field£i2&. Fig.|S{b) reports the raw con- 
ductance curves of the same contact as in Fig.|SJa), mea- 
sured in a magnetic field of increasing intensity. Again, 
the thick line corresponds to the normal-state conduc- 
tance curve obtained here at i? = 6 T (notice that it is 
identical to that of panel (a)). It is clearly seen that, at 
lower fields, the conductance peaks shift towards higher 
energies - a behavior that cannot be explained within a 
single-band model and arises from the stronger suppres- 
sion of the TT-band gap by the magnetic fieldS^Si. The 
curves measured with S = and i3 = 2 T are re- 



FIG. 6: Normalized conductance curves (symbols) of differ- 
ent anomalous point contacts on neutron-irradiated Mg^^B2 
polycrystals at T=4.2 K. The curves are labeled with the An- 
dreev critical temperature but the name of the samples 
and the relevant bulk Tc are also indicated. Lines are the 
best-fit curves given by the two-band BTK model. 



ported, after normalization, in Fig. |Sfc) together with 
the relevant two-band BTK fit. The shift of the con- 
ductance peaks is indicated by the two vertical lines. 
The values of the best-fitting parameters are the fol- 
lowing: = 3.38 meV, F^ = 1.35 meV, = 0.74, 
Ao- = 5.00 meV, F,, = 1.20 meV, Z„ = 0.9 for the zero- 
field curve; A^ = 0.6 meV, F^ = 2.05 meV, = 0.74, 
A^ = 4.85 meV, F^ = 2.75 meV, = 0.9 for the curve 
in magnetic field. The weight was taken equal to 0.8, 
as usual in polycrystalline samples. 

Fig. El clearly shows that, apart from the high value 
of , the anomalous contacts present very regular con- 
ductance curves, with a smooth dependence on magnetic 
field and temperature. Other examples of normalized 
conductance curves of anomalous contacts with rang- 
ing from 35 K down to 22 K are reported in Fig. to- 
gether with the relevant two-band BTK fits. The agree- 
ment between experimental data and fitting curves is re- 
markably good. The gap amplitudes extracted from these 
fits (and from the fit of the curves in other anomalous 
contacts) are reported as a function of the local critical 
temperature in Fig. |7| (black symbols). The behavior 
of the gaps measured in standard contacts is reported for 
comparison (open symbols). It is clear that, while the A^ 
vs. curve is rather similar in standard and anomalous 
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FIG. 7: Gap amplitudes from PCS measurements, as a func- 
tion of the Andreev critical temperature, T^, in anomalous 
contacts (T^ > Tc) (black symbols) compared to the gaps in 
standard contacts (open symbols). Solid lines are only guides 
to the eye, while the straight dashed line indicates the value 
of the gap for a conventional, BCS superconductor. 



contacts, the trend of is fairly different. In anomalous 
contacts, the small gap A^r tends to remain constant or 
slightly increases on decreasing T^, which generally indi- 
cates a substantial increase in the interband scattering^iS. 
Extrapolating the experimental curves to lower critical 
temperatures suggests that the two gaps might tend to a 
common value of about 3 meV and reach it at a critical 
temperature of about 18 K, as in C-doped MgB2 sin- 
gle crystals^. Finally, it is worthwhile to note that, in 
anomalous contacts, the BCS rule for the gap ratio is no 
longer violated. 

The experimental trend of the gaps in anomalous con- 
tacts as a function of can be analyzed within the 
two-band Eliashberg theory, as we did for standard con- 
tacts. Again, the first and simplest possibility consists in 
keeping all the parameters fixed to their values in pris- 
tine MgB2 and simply increasing the interband scatter- 
ing rate Fo-Tr to simulate the disorder due to irradiation. 
The theoretical curves are reported in Fig.lSJa) as dashed 
lines, and the corresponding values of Fo-tt are the same 
we already showed in Fig. 01b). As in standard con- 
tacts, the experimental values of A^ are incompatible 
with this simple picture. As we did for standard con- 
tacts, the next step towards a theoretical reproduction 
of the experimental data is to allow variations in the a- 
band density of states at the Fermi level, neglecting the 
increase in interband scattering. The best fit of the ex- 
perimental gaps versus Tc is obtained by decreasing al- 
most linearly Na-{E-p) from 0.30 down to 0.23 states/(eV 
unit cell) while ranges from 38.8 K to 20 K. The re- 
sulting curves, reported as thin solid lines in panel (a), 
clearly do not follow the experimental values of A^. A 
much better result can be obtained by using both Ncr{Ep) 
and Fq-tt as adjustable parameters to fit the experimen- 
tal data. The best-fitting A„{T^) and A^{T^) curves 
are reported in Fig. |HIa) (thick solid lines) while the val- 
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FIG. 8: (a) Gap amplitudes from PCS measurements in 
anomalous contacts (symbols), as a function of the Andreev 
critical temperature T^. Open and filled symbols refer to 
sample P4 and P6, respectively. Lines: theoretical curves ob- 
tained within the Eliashberg theory by using as adjustable 
parameter only Fo-n- (dashed lines), only NcriEp) (thin solid 
lines) or both of them (thick solid lines). The values of Fo-vr 
and NcriEp) that give the thick solid lines in panel (a) are 
reported as a function of the critical temperature in panels 
(b) and (c), respectively. 



ues of the relevant parameters are reported as a function 
of in Fig. ^h) and (c). The experimental gap val- 
ues measured in anomalous contacts look to be quite well 
reproduced by taking into account the reduction in the a- 
band DOS and the increase in interband scattering. The 
interesting point is that the data in anomalous contacts 
cannot be fitted without interband scattering - while the 
data in standard contacts can, as clearly shown in Fig2| 
This is a consequence of the different behaviour of the 
small gap in the two cases and suggests that standard 
and anomalous contacts occur in regions of the sample 
with different degrees of disorder. 



V. POSSIBLE ORIGIN OF ANOMALOUS 
CONTACTS 

Let us summarize here the properties of anomalous 
contacts that follow from the above. 

1) Their critical temperature is much greater than 
the bulk Tc (from 10 to 25 K more). 

2) In these contacts Ao- is very similar to that measured 
in standard contacts with the same T^, while A^^ is a 
little greater indicating a possible enhancement of inter- 
band scattering. 

3) The probability of finding anomalous contacts in- 
creases on increasing neutron irradiation. 

Properties (1) and (2) might indicate that anomalous 
contacts are established in regions of the sample that 
are either less damaged or partially "reconstructed" , but 
in any case less disordered than the surrounding mate- 
rial. Point (3), however, rules out the former possibil- 
ity so that anomalous contacts are most probably due 
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to reconstructed regions. It has been shown in a recent 
paperiSI that heavily neutron-irradiated MgB2 (with 
as small as 5 K) thermally annealed at sufficiently high 
temperatures (up to 500 °C) and long times (24 h) al- 
most recovers all the characteristics (i.e. cell parameters, 
critical temperature, residual resistivity) of the pristine 
samples. In our case, there is no thermal treatment of 
the samples after irradiation, so one can wonder whether 
similar annealing effects can be due to the PCS measure- 
ments - i.e., because of the current locally injected in the 
sample through the point contacts - or from irradiation 
itself above a certain threshold dose. 

Let us analyze first the hypothesis that local annealing 
occurs due to the PCS technique. To do so, let us focus 
for convenience on the contact made on sample P6 whose 
conductance curves are shown in Fig. ^ The normal- 
state resistance of this contact is i?N = 40 and, as 
already said, the shape of the curves (i.e. the presence 
of small dips^^ and the offset of the curves on increas- 
ing temperature) tells us that the contact must be in 
the diffusive regimoii. If only one contact was estab- 
lished between sample and counterelectrode, its radius 
a - evaluated from the resistance i?N by means of the 
Wexler formula (eq. ^ - would be a ~ 90 A, which has 
to be compared with the effective mean free path ~ 5 A 
and with the coherence length S, — 100 A-'. None of the 
two conditions for spectroscopic analysis to be possible, 
i.e. a £ (ballistic conduction) and a <C would be 
fulfilled. In this situation, at bias voltages comparable 
to the energy gap, the carrier velocity would largely ex- 
ceed the depairing valueS^ and superconductivity would 
be destroyed in a region close to the contact of radius 
three times larger than ^, with consequent loss of the 
Andreev signal. Of course this contrasts with the evi- 
dence of spectroscopic information present in the curves 
of Fig. n (c). Moreover, in these conditions the heating 
in the contact would be greater than experimentally ob- 
served. In fact, using the standard equation valid for a 
circular apertureii 



rTi2 

max bath 



(2) 



one obtains that, for a total voltage drop of the order 
of Vpcak, the temperature of the contact would reach the 
bulk Tc = 8.7 K when Tbath = 6.5 K. This contrasts 
with the curves reported in Fig. ^ where the difference 
between Tc and is less than 0.5 K. 

To reconcile the experimental findings with the value of 
the contact resistance, we are forced to admit that more 
than one contact is established between sample and coun- 
terelectrode. Actually, this is a rather natural assump- 
tion, considering the nature of our point contacts whose 
macroscopic area is about 2000 jixD? . In our case, it can 
be shown that the existence of w 20 diffusive contacts 
(supposed identical and with normal resistance i?N = 800 
VL) can perfectly explain the observed heating in the con- 
tact region. In this case one obtains, for each contact^ 
a=8.2 A, and the temperature in the contact reaches the 
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FIG. 9: Schematic representation of a point contact of radius 
a with the region of radius rj driven normal by the current 
density, which exceeds the critical value in the vicinity of the 
striction. Points A and B are defined in the text. 



bulk Tc for V = Fpcak when Tbath = 8.2 K. Moreover, 
when the bias is of the order of magnitude of the gap, 
the current density j in the contact region is overcritical 
(take into account that jc — 2 ■ 10^ A/cm^ in P6^"'^) but 
the distance rj in the superconductor over which j de- 
cays to the critical value is of the order of 80 A< ^ thus 
ensuring the spectroscopic properties of the contact. 

It turns out from the above that the current injected 
during the measurement of the I — V curve gives rise to 
little or no heating in the contact region, even in the worst 
contact we measured. However, it is worth recalling that, 
whenever necessary, we tuned the normal-state resistance 
by applying voltage pulses of some Volts for ~ 20-^80 ms, 
as experimentally determined. At these bias values, rj is 
certainly greater than ^ so that a big normal region is 
formed (see Fig. ^ in which a very intense current flows 
for a few tens of milliseconds. The normal region and 
the contact itself are then quickly heated above the bath 
temperature (Tbath = 4.2 K). It can be shown that tem- 
peratures of several hundred Kelvin are easily reached 
in the contact. This is witnessed, for example, by the 
early observation of an anomaly at about 250 mV in the 
d'^V/dP of Fe-Fe homocontacts, associated with the fer- 
romagnetic transition of iron at the Curie temperature 
Tc = 770 10^28 

To evaluate the maximum temperature reached in our 
case, let V be the total measured voltage drop between 
sample and counterelectrode and V — Vq the voltage drop 
in the contact itself (let us refer again to the contact on 
sample P6 for convenience), so that Vq is the potential 
difference across the (hemispherical) normal region of ra- 
dius Tj (see Fig. O. If jg is the current density in the 
orifice, then at a distance r from it one has 



jir) = jo 



a 

2^' 



(3) 



The current flowing through the contact is / = (y — 
Vq) /Rpc- Here, Rpc is the resistance of the point contact 
according to eq. ^ where, in the Maxwell term, p{T) is 
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the resistivity of the sample in the normal state, at the 
temperature it will reach in the contact region because of 
the Joule effect. Let us call pavc the average of p{T) in the 
temperature range to be determined (in sample P6 the 
RRR is so small that the value pavo ~ 160 /ii7cm can be 
acceptable for a wide range of temperatures). Using these 
expressions one can calculate the voltage drop across the 
normal region that is given by 



Vo = — where 



\ + K 



K 



1 1 



2ttRpc \a r\ 
Using eq. |3| one also obtains 



l {V-Vo) 
27ri?pcic ' 



(4) 



(5) 



In our case (20 contacts, each with a = 8.2 A) and for 
V =1Y, the solution of Eqs.H and gives Vo = 0.257 V 
and Tj — 2576 A. The maximum temperature in the con- 
tact region (point A in fig. EJ, evaluated from eq. [5]and 
from the voltage drop in the Maxwell part of the contact 
is Ta — 1300 K. The maximum temperature reached in 
the normal region of radius rj can be evaluated by ask- 
ing that the thermal energy generated within the normal 
volume by heating effects equals the flux of heat current 
across the boundary. After suitably simplifying the com- 
plex geometry of the problem, we estimate the maximum 
temperature at the center of the normal region (point B 
in fig. O to be Tb ~ 640 K. An alternative approach to 
evaluate Tb consists in using eq.|21that, even if originally 
derived for a circular aperture^^, approximately holds in 
this case too. This approach gives the result Tb ~ 820 
K. These estimated temperatures are higher than those 
used by Wilke et ali21, even if the duration of the heat- 
ing process is by far shorter. According to our results, 
annealing processes are most probable in the contact re- 
gion, i.e. close to the physical interface between the two 
materials. 

The hypothesis of local annealing as the origin of 
anomalous contacts is thus very reasonable and well 
rooted in the physics of point contact spectroscopy. How- 
ever, we examined as well the other possibility, i.e. pre- 
existent regions with higher Tc than the surrounding ma- 
terial, originated by local reconstruction due to irradia- 
tion. Annealing effects due to the irradiation itself can 
indeed occur, if irradiation takes place at low temper- 
ature, because of the stimulated recombination of close 
Frenkel pairs?^. If, otherwise, irradiation is carried out at 
room temperature or above, competing phenomena such 
as creation and annihilation of point defects or formation 
and coagulation of defect clusters can partly compensate 
each other--, giving rise to saturation in some physical 
parameters. Similar phenomena are suggested, in our 
case, by the tendency to saturation in po: and the 
c-axis parameter™ at very high neutron doses - being 



the other possible reason of saturation, i.e. the complete 
amorphization, ruled out by the sharpness of the X-ray 
peaks. Whatever the exact nature of the reconstruction 
process, locally-annealed regions should feature higher 
critical temperature than the remaining part of the sam- 
ple, but since their presence is neither detected by suscep- 
tibility, nor by specific heat and resistivity measurements, 
they should represent a negligible part of the sample vol- 
ume and should be imagined as isolated regions of small 
size. Moreover, if the macroscopic correlation between pQ 
and Tc observed in irradiated and annealed samplesiiiSi 
is to be conserved also on a local scale, these regions 
are expected to feature higher conductivity and greater 
density of states than the surrounding matrix. 

With these ideas in mind, we performed room- 
temperature SEM and STM analysis of the most irradi- 
ated sample, where the probability of finding anomalous 
contacts was the highest. FESEM morphological images 
of sample P6 showed large, well connected grains with 
smooth surfaces. Microprobe analysis (EDX with SiLi 
detector sensitive to light elements, B included) showed 
no trace of chemical species other than Mg and B. Owing 
to the relatively large size of the grains, we were able to 
perform scanning tunneling microscopy at room temper- 
ature on the surface of a single grain. The topographical 
image reported in Fig. llOf a') shows rather smooth mod- 
ulations on a length scale ~ 10 nm in the xy plane and 
small brighter "dots" that look like protrusions. How- 
ever, morphological SEM images on the surface of grains 
on a similar scale (100 nm x 100 nm) show no trace of 
such spots. The reason is that STM is sensitive not only 
to morphology, but also to the local density of states; re- 
gions with higher DOS appear brighter in the STM map 
because they give rise to a higher conductance across 
the tunnel junction. The topographical signal can be 
removed by operating in STS mode, i.e. keeping the tip- 
to-sample distance constant. The resulting map of the 
current measured across the junction at a constant bias of 
0.1 V is shown, in inverted gray scale, in Fig. lior b). The 
darker spots (corresponding to higher currents) clearly 
correspond to the bright spots of Fig. llOf a). while the 
smooth modulation is no longer observed. For further 
confirmation. Fig. Iior d reports the I-V characteristics 
measured in two points (indicated by small squares in 
panels (a) an (b)) together with the average of the I-V 
curves measured in the whole region. It is clear that the 
"dots" observed in the STM maps represent very small 
regions (typically ~ 1 nm) with higher conductivity 
with respect to the surrounding material, that is exactly 
what one would expect for locally annealed regions with 
higher critical temperature than the bulk. 

The problem now arises of understanding if these dots 
can be superconducting above the critical temperature of 
the matrix in spite of their very small size (if compared to 
the bulk coherence length ^ ~ 10 nm). A large number of 
experimental and theoretical papers have been devoted 
to the so-called "size effect" in isolated superconduct- 
ing nanoparticles2i*22i^^^ but very little is known about 
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FIG. 10: (a) Scanning tunnel microscopy (STM) map of the surface of a single grain in sample P6, measured at room 
temperature. The map was taken with a fixed current equal to 2.0 nA. (b) Map of the current across the N-I-N junction 
(tip/air gap/sample) measured at constant tip-to-sample distance and constant voltage equal to 0.1 V. (c) I-V characteristics 
of the N-I-N tunnel junction in points 1 and 2 of the topographical image (a) . The average 1-V characteristic is also shown for 
comparison. 



the case in which these particles are embedded in a con- 
ductive matrix. To the best of our knowledge, only one 
experimental investigation was reported showing a re- 
duction in the critical temperature of lead nanoparticles 
embedded in a metallic matrix on decreasing the particle 
sizeS. The reduction in becomes effective below 20 
nm (to be compared with the bulk ^ ~ 90 nm) but is 
much smoother than for Pb isolated particles^fe^ so that 
even for very small size (~ 5 nm) is only reduced by 
30% with respect to the bulk value. On this basis, we can 
argue that regions of 3-5 nm in size with a partially recon- 
structed lattice (which would be superconducting even if 
isolated, according to Anderson's criterion^^) could well 
be superconducting even above the critical temperature 
of the matrix and possibly give rise to proximity effect on 
the surrounding normal material. Similar regions, made 
up of clusters of nanoscopic "dots" , are indeed observed 
by STM in some part of the grain surface. Their density 
is much lower than that of the bright dots in Fig. Iior a) 
and this raises the problem of explaining the very high 
probability of anomalous contacts. Again, the answer can 
be given by the technique we used to tune the contact 
characteristics. Owing to their higher conductivity with 
respect to the surrounding matrix, these regions could in 
fact be privileged for the formation of new conduction 
channels when a voltage pulse is applied. The resulting 
new contact would then be dominated by the conductiv- 
ity (and the critical temperature) of these regions. 

Both the proposed mechanisms of formation of anoma- 
lous contacts require the application of a voltage pulse, 
either to provoke the local annealing in the contact region 
or to select the preexistent regions with higher conduc- 
tivity (and possibly higher local Tc). Looking for indirect 
support to our hypotheses, we checked all the contacts 
we studied during several months and we realized that 



indeed only "modified" contacts show anomalous charac- 
teristics and the few standard contacts we were able to 
obtain in highly irradiated samples (among which the one 
shown in Fig.^c)), were actually "as-made". This evi- 
dence further supports our picture and indicate that one 
of the mechanisms described above could really explain 
the origin of the anomalous contacts. 

The local annealing induced by our PCS technique ap- 
pears to be the best understood and most likely pro- 
cess giving rise to anomalous contacts. It easily accounts 
for all the experimental facts, i.e.: i) the occurrence 
of anomalous contacts only after voltage pulses; ii) the 
increasing-with-fluence probability to find such contacts 
(related to the greater concentration of defects that can 
annihilate on annealing); iii) the different behaviour of 
with respect to standard contacts (due to the per- 
sistence of additional disorder in the annealed regions, 
which is consistent with the partial recovery of the pris- 
tine properties under annealing^^). 

The second picture, in which anomalous contacts 
are established on pre-existent regions partially recon- 
structed by irradiation requires making some hypotheses: 
i) the regions with higher DOS observed by STM are less 
disordered than the surrounding matrix; ii) regions of 3-5 
nm in size can develop superconductivity even above the 
bulk Tc, iii) the density of these regions is so small that 
the probability for them to occur in "as-made" contacts 
is vanishingly small; iv) these regions are privileged for 
the formation of new conduction channels when a voltage 
pulse is applied. Some of these points deserve further in- 
vestigation, also for a better understanding of the nature 
of defects in irradiated MgB2. 

Nevertheless, at this stage of investigation, it seems 
very unlikely that the nanoscale inhomogeneities do not 
play any role in the formation of the anomalous con- 
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tacts; on the other hand, the mechanism of local current- 
induced annealing appears very convincing. The most 
reasonable picture is then that these two effects coexist 
and interact. On one hand, local heating in the contact 
region might help the migration of defects and their clus- 
tering together with the partial re-arrangement of nuclei 
in the lattice disordered by irradiation. On the other 
hand, more conductive regions might really be preferred 
channels for the current flow across the junction when 
the voltage pulse is applied, and thus become the centers 
from which the annealing process starts. 

VI. CONCLUSIONS 

In conclusion, we have presented the results of point- 
contact measurements in polycrystalline samples of 
neutron-irradiated Mg-'^'^B2. By using the soft point- 
contact technique developed for MgB2 and related com- 
pounds, we measured the dependence of the gaps on the 
local critical temperature of the contacts, T^. The re- 
sulting trend is in very good agreement with the results of 
specific-heat measurements, especially in the most irradi- 
ated samples, and perfectly confirms the first observation 
of gap merging in undoped MgB2^. This is particularly 
noticeable since the two techniques are completely differ- 
ent and probe the surface and the bulk of the samples, 
respectively. An analysis of the experimental gap trend 
within the Eliashberg theory shows that a major role is 
probably played by the decrease in the density of B pxy 
states, even if an increase in interband scattering can be 
present as theoretically expected. A fit of the gaps is 
however not possible in the whole range of critical tem- 



peratures because in a certain intermediate region both 
the gaps are smaller than the BCS value. 

A striking experimental result was the occurrence of 
anomalous contacts with higher than the bulk Tc^ 
in a percentage increasing with fluence and approaching 
100% in the most irradiated sample. The conductance 
curves of these contacts are perfectly fitted by the two- 
band BTK model, and their temperature and magnetic- 
field dependencies are perfectly standard. However, the 
trend of the small gap A^r extracted from their fit differs 
from that obtained in standard contacts and can be in- 
terpreted within the Eliashberg theory as being due to a 
more effective interband scattering. Annealing effects, ei- 
ther due to our particular PCS technique or to irradiation 
itself, have been proposed to explain this anomaly. The 
two pictures have been carefully investigated, both the- 
oretically and experimentally. The first one - in which 
local annealing arises from the voltage pulses we used 
to tune the contact resistance - appears to be the most 
likely, but the observation by STM of nanoscale regions 
with higher DOS than the surrounding matrix could sup- 
port the second as well. Actually, a cooperative interac- 
tion of the two phenomena looks very probable at this 
stage of investigation. 

We are indebted to Ruggero Vaglio for suggestions 
about the intepretation of anomalous contacts and to 
Pratap Raychaudhuri for enlightening discussion about 
heating in the contact region. Special thanks to Marina 
Putti for fruitful collaboration and help. This work was 
done within the PRIN Project No. 2004022024 and the 
INTAS Project No. 01-0617. V.A.S. acknowledges sup- 
port by Russian Foundation for Basic Research (Proj. 
No. 06-02-16490). 
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